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Improvement  of  Real  Time  foFj 
Prediction  for  OTH  Radars 


1.  INTRODUCTION 

Dandekar  and  Buchau^  (1986)  proposed  a  technique  for  improving  predictions  over  an 
interval  of  6  hours  around  the  sunrise  transition  period.  The  technique  uses  the  hourly  average 
gradient  from  the  previous  four  days’  of  foFj  observations  and  the  currently  (present)  observed  foF2 
to  predict  f^Fj  for  the  next  hour.  An  additional  improvement  is  achieved  by  adding  a  second  order 
correction  to  the  gradient  by  considering  the  sign  of  the  change  (  positive/negative )  of  the  present 
f^jFj  with  respect  to  the  prior  four  days’  average  f^Fj  for  the  same  hour.  In  this  report  the  same 
technique  is  tested  for  24-hour  real  time  prediction  of  fc)F2'  Thus,  with  a  minor  modification,  the  old 
scheme  is  found  to  be  very  promising. 

For  a  simulation  of  the  OTH  East  Coast  radar  situation  three  high  latitude  ionospheric  stations 
in  Canada:  St.  Johns  (Newfoundland),  Ottawa,  and  Winnipeg,  are  used.  To  cover  different  levels  of 
solar  activity,  foF2  data  are  selected  for  high  (calendar  year  1969)  and  low  (calendar  year  1976)  solar 
activity.  The  use  of  the  algorithm  proposed  herein  improves  the  hourly  prediction  by  45  percent  over 
the  predictions  obtained  from  using  just  the  lONCAP  modeP  (Lloyd  et  al,  1978). 

The  first  section  presents  the  algorithm  and  the  analysis  procedure  used  in  the  study.  Section 
2  discusses  the  distribution  seen  in  the  foF2  data  base,  which  forms  the  basis  for  the  algorithm. 
Section  3  presents  results  using  the  algorithm.  The  last  section  summarizes  the  merits  of  the 
proposed  scheme  for  improving  foF2  predictions. 
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1.1  Algorithm  for  the  Prediction  of  foFj 

The  proposed  algorithm  is  given  by  Eq.  (1), 


is  (4*’: 
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where  the  term  on  the  left  hand  side  of  Eq.  (1)  is  the  predicted  for  the  n+1  st  day  (today)  for  the 
next  hour,  T  is  the  'present'  hour,  and  n  refers  to  the  prior  'n'  days.  The  first  term  on  the  right  hand 
side  (RHS)  is  the  observed  foF2  today  (n+1)  for  the  present  hour(i). 

The  second  term  (with  the  summation  sign)  is  the  averaged  change  in  foFj  between  the  i+1 
and  'i'th  hour  fi-om  the  prior  n  (  n=4,  used  in  this  study )  days.  The  term  f^Fj  (with  overline)  is  the  fJFj 
averaged  for  n  prior  days  for  the  'i'th  hour  and  is  given  by  Eq.  (2).  Also,  it  can  be  used  as  the 
prediction  for  the  present  ('i'th)  hour  for  (the  n+lst  day)  today,  if  the  'present'  foFj  (the  first  term  on 
the  RHS)  is  not  avdlable.  In  Eq.  (1)  the  last  term  in  the  end  bracket  is  a  second  order  correction  to 
the  slope  where  a  positive  or  a  negative  sign  is  used,  depending  on  the  empirically  determined  time 

intervals. 

Thus  the  prediction  of  foFj  consists  of  the  'present'  (  observed)  foFj,  plus  the  change  in  f^Fj 
over  the  next  hour  averaged  for  prior  'n'  days,  and  with  a  correction  term  that  modifies  the  slope  of 

foFa- 

The  estimate  of  today's  can  be  obtained  from  the  average  of  previous  n  days'  f^F^  and 
is  given  by  Eq.  (2), 


(2) 


In  place  of  the  observed,  present  (day  n+1,  'i'th  hour)  foF2,  one  can  compute  this  first 
term  on  the  RHS  of  Eq.  (1),  as  the  averaged  f/2  from  Eq.  (2).  A  similar  form  of  weighted  mean  was 
proposed  by  Rush  and  Gibbs'  for  f„F2  predictions  in  1973.  In  their  scheme  the  prior  day  foF2  was 
multiplied  by  a  weight  factor  that  decreased  with  the  order  of  the  day  (  prior  day  1;  factor  n,  prior  day 
2;  factor  n-1,  prior  day  3;  factor  n-3,..  up  to  5  days ).  In  our  scheme  (Eq.  (1))  we  compute  a  simple 
average  of  foF2  without  any  such  weight  factors. 
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The  proposed  method  is  shown  schematically  in  Figure  1  for  a  hypothetical  set  of  data.  The 
figure  is  divided  into  two  sections.  The  left  hand  section  (LHS)  is  for  the  period  of  increasing  fj^2 
and  the  right  hand  section  (RHS)  is  for  the  period  of  decreasing  f^Fj- 1*'  two  sections  of  the 
figure,  we  see  that  for  a  given  trend  (increasing  or  decreasing  f^j)  the  change  of  (positive  to 
negative)  sign  of  the  correction  term  in  the  last  bracket  on  the  RHS  of  Eq.  (1)  produces  an  opposite 
(converging  or  diverging)  change.  That  is,  for  the  next  hour  deviates  less  (converges)  or  more 
(diverges)  fi'om  the  first  order  prediction.  The  comparison  of  the  two  periods  shows  that  the  same 
sign  produces  an  opposite  change  in  these  trends.  Also  note  that  in  Figure  1  on  the  last  line  of  the 
algorithm,  the  numerator  is  a  difference  term.  An  interchange  in  the  terms  in  the  numerator 
interchanges  the  role  of  the  positive  and  negative  sign  for  this  term.  In  the  figure  the  lines  closer  to 
the  central  line  show  a  conver^g  trend  and  the  lines  farthest  from  it  indicate  a  diverging  trend.  Thus 
we  would  like  to  know  which  trend,  converging  or  diverging,  produces  better  hourly  predictions  of 
f^2>  ^  given  time. 

We  propose  (later  in  the  discussion)  the  use  of  the  scheme  in  the  left  hand  side  of  Figure  1 
for  04-16  LT  and  the  right  hand  side  of  Figure  1  for  16-04  LT  (converging  trend  in  both  cases)  for 
predictions  of  foF2. 

1.2  Other  Schemes 

The  lONCAP  method  is  used  as  the  baseline  for  determining  the  improvement  achieved  by 
any  of  these  methods.  Two  schemes  that  truncate  the  RHS  of  Eq.  (1)  at  different  terms,  and  a 
correction  term  of  opposite  sign  in  the  modified  slope  method,  were  also  tested.  The  use  of  the 
equivalent  of  a  perfect  prediction  from  lONCAP  allowed  a  test  of  how  well  a  perfect  prediction  from 
lONCAP  would  serve  an  operational  OTH  radar.  It  was  found  that  even  a  perfect  prediction  (of  the 
monthly  median  f^Fj  for  each  hour)  from  lONCAP  does  not  improve  the  short  term  prediction  of 
foF2  very  much. 

In  the  analysis  the  mean  monthly  sunspot  number  (published  by  National  Oceanic  and 
Atmospheric  Administration  (NOAA),  Boulder,  Colorado,  US)  is  used  with  lONCAP.  In  the  OTH 
operation.  Air  Weather  Service  (AWS)  computes  the  effective  sunspot  number  from  foF2  observed 
from  the  network  of  50  globally  spaced  ionosonde  stations.  The  AWS  scheme  of  the  effective 
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decreasing  trend  for  these  sections. 


sunspot  number  could  not  be  used  in  this  analysis.  Note  that  as  a  starting  point,  the  AWS  uses  the 
ITS-78  model  (Barghausen  et  al*  1969  ),  which  is  the  precursor  of  the  lONCAP  model. 

The  monthly  median  predictions  are  computed  for  each  hour  from  lONCAP  using  the 
observed  mean  sunspot  number.  These  predictions  are  then  used  to  determine  the  improvement  in 
prediction  achieved  by  the  schemes  described  above.  Thus  we  compare  results  from  six  methods. 

1)  lONCAP, 

2)  observed  medians  /accurate  lONCAP, 

3)  prior  four  days'  averaged  foFj  for  the  next  hour, 

4)  the  averaged  slope  method, 

5)  modified  slope  method  using  converging  trend,  and 

6)  modified  slope  method  using  diverging  trend  (depending  on  the  sign  used  in  the  last  bracket 
ofEq.  (1)). 

Three  examples  are  presented  to  demonstrate  the  improvement  of  the  predictions  by  use  of 
the  algorithm.  Figure  2  shows  a  mass  plot  of  foF2  data  from  Argentia,  Canada  for  November  1990 
along  with  the  observed  median  values  and  the  lONCAP  predictions  (SSN=140).  Note  that  the 
observed  medians  nicely  follow  the  lONCAP  predictions,  but  the  lONCAP  predictions  on  the 
average  are  1  MHz  greater  than  the  median  values.  Also,  the  median/IONCAP  provides  a  single 
value  for  f^Fj  for  each  hour  of  a  given  month,  whereas  the  observations  cover  a  wide  range  of  2.5 
MHz  at  09  UT,  to  14.2  MHz  at  17  UT.  Note  that  the  observations  for  Julian  day  90-331  (November 
27,1990,  with  a  major  magnetic  storm)  show  a  large  difference  from  the  remaining  data.  Figure  3 
shows  the  predictions  for  Argentia  from  the  full  algorithm  of  Eq.  (1)  along  with  the  observations 
for  three  days  selected  from  Figure  2.  Figure  3  shows  that  the  predictions  from  the  algorithm  are 
close  to  the  actual  observations  and  show  an  improvement  in  the  prediction  over  both  the  observed 
median  and  lONCAP  predictions.  Considering  the  large  deviation  of  the  frequency  for  10-24  UT  on 
Julian  day  33 1  the  predictions  from  the  algorithm  are  very  good.  This  demonstrated  success  of  the 
algorithm  must  be  validated  with  a  larger  data  base. 

Before  applying  the  algorithm  to  a  larger  database,  the  systematic  behavior  of  the  f^Fj 
observations  is  examined.  This  systematic  behavior  allows  the  use  of  the  proposed  algorithm  for  a 
significant  reduction  in  the  prediction  error. 
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Figure  2.  Mass  plot  of 
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Figure  3.  Comparison  of  observations  with  predictions  from  the  proposed  algorithm. 


2.  DATA  BASE 

The  ionospheric  stations  representing  the  northern  sector  of  the  OTH  East  Coast  Radar  are 
listed  in  Table  1.  The  selection  of  the  stations  is  very  much  dependent  on  the  availability  of  a 
reasonably  continuous  set  of  observations  of  f<,F2  over  a  long  time  period.  The  table  lists  the  name, 
geographic,  and  corrected  geomagnetic  (CG)  coordinates  of  the  stations.  St.  Johns,  Canada  is  in  the 
northern  sector  of  the  OTH  ECRS  coverage.  Most  often  this  station  is  south  of  the  midlatitude  F- 
layer  trough.  Although  the  second  station,  Ottawa,  Canada,  does  not  lie  in  the  OTH  coverage,  its 
geograpWc  latitude  is  a  good  representation  of  the  OTH  coverage  area.  The  third  station,  Winnipeg, 
Canada,  is  both  geographically  and  geomagnetically  north  of  both  Ottawa  and  St.  Johns,  and  is  close 
to  the  auroral  oval.  The  ionospheric  characteristics  of  Winnipeg  are  close  to  those  of  Goose  Bay, 
Canada.  Goose  Bay  lies  in  the  northern  sector  of  the  ECRS.  Because  Goose  Bay  did  not  operate 
prior  to  1972  it  could  not  be  used  in  this  study  to  cover  a  full  solar  cycle. 


Table  1.  Ionospheric  Stations  Used  in  the  Study 


Geographic 

Corr.  Geomagnetic 

Station 

Lat.N 

Long.  E 

Lat.N 

Long.  E 

St.  Johns 

47.6 

307.3 

57.6 

29.1 

Ottawa 

45.1 

283.8 

58.5 

355.7 

Winnipeg 

49.8 

265.6 

61.1 

323.5 

For  the  selected  stations,  the  foFj  data  for  the  calendar  years  1969  and  1976  are  used  as 
representative  of  the  high  and  low  phases  of  solar  activity.  The  observed  averaged  monthly  sunspot 
numbers  (SSN)  for  the  years  of  high  and  low  solar  activity  are  106  and  1 1  respectively.  Examination 
of  the  diurnal  behavior  of  f^Fj  shows  why  the  proposed  algorithm  provides  good  fJ^j  predictions. 
This  is  shown  in  Figures  4A-E  for  1969  (no  continuous  data  are  available  for  Winnipeg)  and  1976. 
Each  figure  is  divided  into  12  blocks.  Each  block  presents  the  month  (noted  in  the  upper  comer  of 
the  block).  For  each  month  there  are  three  curves  showing  the  upper  quartile,  the  median,  and  the 
lower  quartile.  Gaps  in  the  data  are  shown  by  dotted  lines.  The  points  to  note  are:  a)  the  diurnal 
variation  of  f^Fj  shows  a  maximum  around  1500  LT  and  a  minimum  in  the  predawn  hours. 
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Figure  4A.  Diurnal  variation  of  at  St.  Johns  in  1969. 


Figure  4B.  Diurnal  variation  of  at  Ottawa  in  1969. 
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Figure  4C.  Diurnal  variation  of  at  St.  Johns  in  1976. 
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Figure  4D.  Diurnal  variation  of  at  Ottawa  in  1976. 
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Figure  4E.  Diurnal  variation  of  at  Winnipeg  in  1976. 


b)  between  the  maximum  and  the  minimum  the  diurnal  variation  changes  by  a  factor  of  2,  c)  the 
diurnal  variation  is  largest  in  winter  months  (October  to  March)  as  shown  in  rows  1  and  4  and 
relatively  smaller  in  summer  months  (April  to  September)  as  shown  in  rows  2  and  3,  and  d)  the  peak 
fjjFj  decreases  by  a  factor  of  2  between  the  high  (1969)  and  low  (1976)  solar  activity  periods.  This 
systematic  diurnal  variation  is  the  reason  that  good  predictions  using  the  algorithm  in  Eq.  (1)  are 
possible. 

The  percent  change  (slope)  in  hourly  f^Fj  is  computed  from  Eq.  (3) 

p  =100 - - i— . (3) 


The  median  hourly  slopes  in  f^Fj  are  shown  in  Figures  5A-E.  Each  figure  presents  the 
contours  of  median  hourly  slopes  for  the  whole  year  for  each  month  for  all  24  hours.  The  positive 
sign  refers  to  an  increase  in  f^Fj  and  the  negative  sign  refers  to  a  decrease  in  f^Fj.  There  are  two 
periods  of  large  changes  (^20  percent),  one  around  the  sunrise  and  the  other  around  the  sunset.  The 
sharpest  increase  in  f^Fj  is  seen  between  04-08  LT.  The  sunrise  changes  are  the  largest  (^40  percent) 
in  the  winter  months. 

The  other  factor  to  consider  is  the  spread  of  f^Fj  for  a  given  LT  hour.  For  each  hour  of  a 
given  month  the  normalized  spread  (d/Z’^)  is  computed  from  the  equation 

fF.  -fF. 

5/ F  =100 - ^ ^ . (4) 

°  2  2fF, 


The  terms  in  Eq.  (4)  are  median  f(,F2  (not  mean  foF2),  and  foF2  ±o  levels.  If  the  foF2 
distribution  is  normal,  the  median  and  mean  foF2  are  equal.  Also,  the  magnitudes  of  ±o  levels  will 
be  the  same.  In  this  analysis  the  median  (50th  percentile)  and  ±a  (84th  and  16th  percentiles,  that  is, 
68  percent  population  around  the  median)  levels  were  determined  from  the  observed  distribution  of 

foF^. 

The  time  dependence  of  (relative)  spread  in  f^Fj  is  shown  in  Figures  6A-E.  Figures  6A-E  have 
the  same  format  as  Figures  4A-E.  These  figures  show  that  the  spread  of  foF2  is  smaller  for  the 
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Figure  5C  Median  monthly  variation(%)  in  f„F2  at  St.  Johns  in  1976. 
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Figure  5E  Median  monthly  variation(%)  in  foF2  at  Winnipeg  in  1976. 
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L.T. 

Figure  6 A.  Percent  variation  between  ±o  levels  of  at  St.  Johns  in  1969. 
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Figure  6B.  Percent  variation  between  ±o  levels  of  at  Ottawa  in  1969. 
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Figure  6E.  Percent  variation  between  ±o  levels  of  at  Winnipeg  in  1976. 


daytime  hours  and  larger  for  nighttime  and  for  the  transition  periods  of  sunrise  and  sunset.  We  find 
that  the  success  of  the  f^Fj  prediction  algorithm  is  inversely  related  to  this  spread;  that  is,  smaller  the 
spread,  better  is  the  prediction  of  f^Fj.  Note  that  we  are  considering  the  relative  spread  in  f^Fj  (not 
the  absolute  values  of  the  spread).  In  the  following,  the  algorithm  is  used  to  compute  the  predicted 
fj,F2  by  all  the  six  methods  mentioned  above. 

3.  ANALYSIS 

The  predictions  from  the  lONCAP  program  are  used  as  the  reference  for  determining  the 
improvement  in  the  prediction  of  foF2  by  the  algorithm  (Eq.  (1)).  Therefore,  the  first  step  is  to 
determine  the  error  between  the  observed  ^F2  and  the  ^F2  predicted  by  the  lONCAP  program.  The 
monthly  averaged  sunspot  number  listed  in  Table  2  is  used  as  an  input  to  the  lONCAP  program. 


Table  2.  Zurich  Sunspot  Number 


YEAR 

1969 

1976 

MONTH 

JANUARY 

104 

8 

FEBRUARY 

121 

4 

MARCH 

136 

22 

APRIL 

107 

19 

MAY 

120 

12 

JUNE 

106 

12 

JULY 

97 

2 

AUGUST 

98 

16 

SEPTEMBER 

91 

14 

OCTOBER 

96 

21 

NOVEMBER 

94 

5 

DECEMBER 

98 

15 

AVERAGE 

106 

11 

The  comparison  of  the  observed  monthly  median  foF2  with  monthly  median  prediction  of  f^Fj  from 
the  lONCAP  is  shown  in  Figures  7A-E  as  a  function  of  time  of  the  day  and  season.  These  figures 
show  that  at  high  solar  activity  (year  1969)  the  range  of  foFj  is  3  to  1 1  MHz,  and  the  standard  error 
between  the  prediction  and  observation  is  less  than  1  MHz,  At  low  sunspot  activity  (year  1976)  the 
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Figure  7A  Monthly  median  contours  of  f^F; 
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Figure 
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at  low  solar  activity  (1976)  at  Ottawa. 
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^2  values  range  from  2  to  6  MHz.  Again  the  error  is  less  than  1  MHz.  At  low  solar  activity 
is  observed  to  be  ^2  MHz  in  the  early  morning  hours  for  much  longer  periods  than  those  predicted 
by  the  lONCAP  model.  As  the  error  at  the  median  level  of  foFj  is  of  the  order  of  1  MHz,  the 
lONCAP  program  cannot  provide  a  real  time  prediction  with  error  <0.5  MHz  sought  by  real  time 
data  collection  systems  like  the  Digital  Ionospheric  Sounding  Systems*  (DISS)  deployed  by  AWS. 

Let  us  consider  the  algorithm  in  Eq.  (1)  for  the  short  term  (next  hour)  prediction  of  foFj.  In 
this  analysis  the  error  in  the  prediction  of  f^Fj  is  computed  from 


AfF=\  fF. 


2 

^prmdiefd  ^ 


(5) 


The  percent  error  (ERR)  in  the  prediction  is  computed  by  using 


ERB{%)  =100 . . (6) 

foF 

^observed 


An  improvement  (IMPR)  in  the  prediction  error  is  computed  by  using  the  error  in  the 
prediction  of  fJFz  the  lONCAP  program  as  the  reference.  It  is  given  by 

JAfPR%)=10C|l - . (7) 

Af  F, 


The  foF2  data  are  analyzed  for  each  hour  on  a  monthly  basis.  The  foF2,  the  percent  error  in 
^2>  the  improvement  I(percent)  are  computed  from  Eqs.  (5)-(7)  for  each  of  the  six  prediction 
methods  listed  in  the  algorithm  section.  From  these,  an  average  value  for  the  respective  error  is 
obtained  for  the  given  hour  of  each  month.  These  values  are  averaged  over  all  hours  to  obtain  a 
monthly  average  error.  Also,  the  values  for  a  given  hour  are  averaged  over  the  year  to  obtmn  an 
hourly  average  error  for  the  year.  All  these  values  are  compared  for  all  the  methods  listed  above  to 
determine  the  merit  of  each  method. 
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The  hourly  fJFj  data  for  the  stations  listed  in  Table  1  are  treated  in  three  categories.  The  first 
category  is  to  use  aU  available  data.  The  second  category  is  based  on  the  observation  accuracy  of  the 
real  time  collection  systems  such  as  the  Automatic  Real  Time  Ionospheric  True  (ARTIST)  height 
profiles  (Reinisch  et  al®  1983)  providing  electron  density  distribution  and  the  layer  parameters.  These 
systems  provide  on-line  foFj  analysis  with  an  accuracy/error  <0.5  MHz  for  fJFj.  Therefore  for  the 
second  category  the  criterion  for  the  error  f^2  k  0.5  MHz  is  used.  The  third  category  is  based  on 
the  ^  of  the  radar  operation  (see  in  the  follovwng  in  the  section  on  relevance  to  the  OTH 
operation).  The  aim  of  the  OTH  radar  is  to  maintain  a  barrier  of  a  given  width  at  a  given  distance 
from  the  radar  at  all  times.  This  is  achieved  through  frequency  management.  For  the  OTH  radar 
every  6  percent  change  in  f^j  sUdes  the  barrier  by  250  km.  Therefore  a  (6*3=)  20  percent  change 
(inability  to  maintain  barrier  of  500  km  if  it  slides  three  times  its  half  width)  in  the  fJFj  is  used  as  the 
error  criterion  in  the  this  category. 

The  data  available  fi-om  the  stations  for  the  error  analysis  are  summarized  in  Table  3.  In  this 
table  the  name  of  the  station  and  the  year  and  the  category  of  f^Fj  data  are  listed  in  the  first  three 
columns.  The  next  24  columns  list  the  number  of  data  points  available  for  each  of  the  three 
categories  specified  above.  The  last  column  lists  the  total  number  of  data  points  available  for  the 
year.  For  each  station  the  first  line  lists  the  number  of  hourly  observations  available  for  the  station. 
The  second  line  provides  the  percent  (number  of  observations  available/365)  hourly  coverage  for  each 
hour.  The  next  three  lines  list  the  percent  amount  of  data  (of  line  1)  used  in  the  respective  categories 
mentioned  on  the  left  hand  side  of  the  table.  The  table  shows  that  the  data  collection  rate  was  high, 
at  about  90  percent,  during  the  high  solar  activity  year  (1969)  and  dropped  down  to  68  percent 
during  the  low  solar  activity  year  (1976)  at  all  three  stations.  The  averaging  scheme  of  the  algorithm 
in  Eq.  (1)  starts  on  the  fifth  day  (n=4)  for  computing  errors  in  f^Fj  by  various  methods.  This  leads 
to  a  reduction  of  13  percent  of  the  useful  data  collected  at  the  stations.  Thus  ideally  in  Table  3  the 
third  line  should  be  87  percent  of  the  first  line.  The  reduced  number  of  predictions :  75  percent 
during  high  solar  activity  period  and  63  percent  during  low  solar  activity  period(  <87  percent)  is  due 
to  discontinuities  (gaps)  in  the  data.  The  next  line  presents  the  percent  data  that  resulted  in  prediction 
errors  in  ^0-5  MHz,  fi-om  at  least  one  of  the  six  methods  with  respect  to  that  from  the  lONCAP. 
The  table  shows  that  at  high  solar  activity,  errors  in  f^Fj  ^0.5  MHz  occur  for  the  lONCAP  predictions 
for  67  percent  of  the  time.  This  occurrence  drops  to  about  45  percent  at  low  solar  activity. 


32 


Table  3.  Number  of  Observations  Available  for  the  Hour 
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The  last  line  of  each  group  shows  that  large  relative  errors  (^20  percent)  in  predictions 
are  least  frequent  (15  percent)  in  the  afternoon.  Their  occurrence  rate  increases  to  50  percent  around 
the  post  midnight  and  early  morning  hours. 

During  the  low  solar  activity  period  of  1976  the  data  collection  at  the  stations  dropped  to  35 
percent  during  the  early  morning  hours.  This  is  also  seen  in  Figures  4C-E  where  the  dotted  line 
segments  indicate  missing  data.  These  periods  form  a  poor  sample  for  the  analysis  of  prediction 
errors  in  f^Fj. 

The  analysis  showed  that  the  modified  slope  method  (Figure  1)  provides  a  slightly  greater 
improvement  (Eq.  (7))  over  the  averaged  slope  method  (algorithm  in  Eq.  (1)  without  the  slope 
correction  term  in  the  square  bracket),  when  the  converging  trend  is  used  for  both  04-16  LT  and  for 
16-04  LT  by  using  the  opposite  signs  for  respective  intervals  for  f^Fj  predictions  from  the  full 
algorithm  ofEq.  (1).  The  results  for  these  two  time  intervals  are  presented  in  Table  4.  In  this  table 
the  first  three  columns  list  the  name  of  the  station,  the  year,  and  the  category  of  the  data  selection. 
The  rest  of  the  table  is  divided  into  two  sections.  In  each  of  these  sections,  there  are  12  columns,  one 
for  each  month,  followed  by  the  column  presenting  the  average  value  for  the  year.  The  aggregate 
improvement  of  the  respective  modified  slope  method  (using  the  corresponding  data  segments)  is 
Osted  in  the  table.  The  division  in  time  intervals  is  not  perfect,  but  on  the  whole  there  is  a  consistent 
improvement  of  3  percent  over  the  averaged  slope  method.  In  Figures  4A-E,  the  change  from  the 
modified  negative  slope  to  the  positive  slope  corresponds  to  the  peak  in  the  diurnal  variation  of  foFj. 
A  look  at  Figure  1  shows  that  the  second  order  correction  using  a  consistent  converging  trend  for 

all  24  hours  results  in  improved  predictions  of  f^Fj. 

The  average  error  (Eqs.  (5)  and  (6))  in  the  prediction  of  foF2  for  each  year  is  computed  over 

the  two  proposed  intervals.  The  results  are  presented  in  Table  5.  The  first  three  columns  on  the  left 
hand  side  list  the  station,  the  year,  and  the  category  for  the  analysis.  The  table  is  divided  into  two 
sections  along  the  columns  for  the  two  time  intervals  specified  above.  In  each  section  the  errors  are 
listed  for  the  sk  prediction  methods.  The  number  of  data  points  and  their  percent  contribution  to  the 
respective  time  intervals  are  also  listed.  The  table  shows  that  the  errors  are  largest  for  the  lONCAP 
method  and  reduce  successively  for  the  other  methods:  observed  median  (which  is  the  accurate 
lONCAP  prediction),  prior  four  days’  average  foF2,  average  slope  method,  and  the  modified  slope 
method.  Note  that  for  one  time  interval  the  positive  slope  method  is  better  and  for  the  other 
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Methods;  1)  loncap,  2)  Observed  Median/  Accurate  loncap,  3)  Prior  Four  Days'  Average,  4)  Average  Slope,  5)  2nd  order  correction 
with  Negative  Sign,  and  6)  2nd  Order  Correction  with  Positive  Sign. 


interval  the  negative  slope  method  is  better,  (essentially  the  improvement  is  associated  vvith  the 
converging  trend  approach).  The  significant  result  is  that  the  error  of  0.9  MHz  from  the  lONCAP 
during  high  solar  activity  is  reduced  to  0.4  MHz  by  the  use  of  the  proposed  scheme.  At  low  solar 
actiwty  the  respective  reduction  is  from  0.5  MHz  to  0.3  MHz.  Thus  the  proposed  method  reduces 
the  error  to  50  percent  of  that  from  the  lONCAP  method.  Considering  the  magnitude  of  errors  (in 
MHz),  significant  improvement  is  observed  at  both  high  as  well  as  low  solar  activity  periods. 

In  terms  of  percent  error,  (^20  percent)  the  large  relative  errors  occur  more  frequently  (65 
percent  of  the  time)  in  the  16-04  LT  interval  than  (35  percent)  in  the  04-16  LT  interval. 

The  percent  improvement  (Eq.  (7))  averaged  over  the  whole  year  is  presented  in  Table  6, 
which  has  the  same  format  as  Table  5.  The  errors  are  averaged  separately  for  the  high  and  low  solar 
activity  periods.  The  improvement  is  better  at  high  solar  activity  (52  percent)  than  at  low  solar 
activity  (35  percent).  Note  that  the  improvements  by  the  observed  median  (that  is,  accurate 
lONCAP)  and  by  'prior  four  days'  average'  methods  are  small;  10  percent  and  14  percent  respectively. 
Thus  a  prediction  fi-om  a)  the  lONCAP,  b)  an  accurate  lONCAP  and  c)  fi-om  the  averaging  method 
(amilar  to  that  proposed  by  Rush  and  Gibbs^  (1973))  are  unable  to  provide  the  reasonably  accurate 
prediction  of  f^j  needed  for  the  real  time  application.  The  division  of  the  algorithm  into  two  time 
intervals  is  slightly  better  than  the  use  of  the  averaged  slope  method  for  improving  the  real  time 
predictions  of  foFj. 

The  averaged  monthly  values  for  the  improvement  using  the  modified  slope  method  for  the 
three  categories  for  the  stations  are  presented  in  Table  7.  The  table  shows  similar  improvements  for 
each  month  in  all  the  categories,  showing  the  improvement  in  the  predictions  based  on  the  algorithm 
inEq.  (1).  The  improvements  in  f^Fj  predictions  are  larger  (52  percent)  at  high  solar  activity  and  a 
little  lower  (33  percent)  at  low  solar  activity.  Ottawa  and  Winnipeg  stations  show  smaller 
improvements  (<25  percent)  in  the  months  of  June,  July,  November  and  December  during  the  period 
of  low  solar  activity. 

The  percent  improvement  is  averaged  for  each  LT  hour  of  the  year  to  check  if  the 
improvement  of  the  foFj  prediction  from  the  algorithm  has  any  time  dependency.  The  results  are 
summarized  in  Table  8.  Figures  5A-E  show  that  the  transition  periods  of  large  slopes  cover  the 
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Table  6.The  Improvement  (percent)  in  the  Prediction  of 
By  Various  Schemes  Compared  to  That  ifrom  lONCAP 


Percent  Improvement 


1 

04-15  LT 

16-03  LT 

METHODS 

METHODS 

STATION 

YR  f„F, 

2 

3 

4 

5 

6 

2 

3 

4 

5 

6 

STJOHNS 

69  ALL 

12 

21 

53 

52 

54 

12 

16 

44 

46 

42 

>.5MHZ 

12 

21 

54 

52 

55 

12 

15 

45 

46 

42 

e>20% 

4 

11 

55 

54 

55 

3 

10 

44 

46 

41 

OTTAWA 

69  ALL 

9 

18 

54 

52 

54 

17 

19 

56 

58 

54 

>.5MHZ 

9 

18 

55 

53 

55 

16 

19 

57 

59 

55 

e>20% 

0 

9 

56 

56 

57 

10 

17 

57 

59 

54 

69  AVG 

8 

16 

55 

53 

55 

12 

16 

51 

52 

48 

STJOHNS 

76  ALL 

11 

11 

27 

25 

29 

17 

14 

32 

35 

28 

>.5MHZ 

10 

11 

26 

24 

29 

16 

12 

31 

35 

27 

e>20% 

8 

10 

28 

27 

31 

20 

16 

34 

38 

30 

OTTAWA 

76  ALL 

13 

15 

36 

35 

38 

17 

16 

39 

42 

36 

>.5MHZ 

11 

14 

34 

33 

37 

16 

14 

39 

43 

37 

e>20% 

8 

19 

36 

37 

40 

16 

15 

40 

44 

37 

WINNIPE( 

j76  all 

12 

15 

38 

32 

39 

6 

4 

21 

24 

18 

>.5MHZ 

9 

14 

38 

32 

39 

4 

2 

18 

21 

14 

e>20% 

9 

18 

40 

37 

43 

-1 

1 

12 

17 

9 

76  AVG 

10 

14 

34 

31 

36 

12 

10 

30 

23 

26 

GR 

AND  AVG 

9 

15 

44 

42 

46 

12 

13 

40 

43 

37 

Methods:  1)  loncap,  2)  Observed  Median/  Accurate  loncap,  3)  Prior  Four  Days'  Average,  4)  Average 
Slope,  5)  2nd  Order  Correction  with  Negative  Sign,  6)  2nd  Order  Correction  with  Positive  Sign. 
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Table  7.  Percent  Improvement  by  Modified  Slope  Method 


month 

STATION 

01  02  03  04  05  06  07  08  09  10  11  12 

AVG 

STJOHNS 

69  ALL 
>.5MHZ 
e>20% 

39  51  52  56  52  44  42  48  64  50  42  49 

39  52  52  56  52  44  44  51  64  51  43  49 

29  61  59  63  61  33  42  43  63  43  47  44 

OTTAWA 

69  ALL 
>.5MHZ 
e>20% 

35  53  56  59  54  64  59  65  65  63  50  40 

35  53  56  60  56  64  60  66  65  64  50  41 

24  63  60  64  58  63  60  66  65  70  48  31 

69  AVG 

34  56  56  60  56  52  51  57  64  57  47  42 

52 

STJOHNS 

76  ALL 
>.5MHZ 
e>20% 

27  29  38  41  33  12  26  25  39  47  21  17 
25  28  38  42  35  11  26  25  39  48  18  12 
25  27  46  42  38  13  21  28  42  55  13  19 

OTTAWA 

76  ALL 
>.5MHZ 
e>20% 

28  27  58  51  41  28  26  29  40  46  30  42 

26  28  58  54  42  26  24  21  42  48  28  43 

33  31  61  61  51  18  10  24  43  54  24  49 

WINNIPEC 

r76  ALL 
>.5MHZ 
e>20% 

29  42  43  34  33  27  18  29  32  43  30  14 

29  44  46  31  32  28  17  30  32  45  29  4 

31  48  52  28  33  26  10  19  32  45  30  13 

76  AVG 

28  34  49  43  38  21  20  26  38  48  25  24 

33 

GRAND  AVG 

31  45  52  51  47  37  36  41  51  52  36  33 
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Table  8.  Annual  Improvement  for  Each  Hour  (Averaged  for  12  Months) 
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periods  04-09  LT  and  16-23  LT.  For  high  solar  activity  (year  1969)  the  improvement  in 
predictions  is  good  for  all  the  LT  hours.  At  low  solar  activity  the  improvement  in  fJFj  prediction  is 
low  (<25  percent).  A  factor  for  such  poor  performance  would  be  the  large  spread  in  fJFj  at  the  given 
hour  seen  in  Figures  6  A-E. 

The  number  of  occurrences  of  poor  performance  (improvement  <25  percent)  for  f^Fj 
predictions  from  the  algorithm  were  counted  for  the  hourly  values  of  the  improvement  (Eq.  (7)).  At 
high  solar  activity  the  algorithm  performs  poorly  for  10  percent  of  the  time,  with  the  occurrence  more 
frequent  (2:1)  in  the  16-04  LT  interval  than  in  the  04-16  LT  interval.  At  low  solar  activity  (1976) 
the  poor  performance  occurs  for  25  percent  of  the  time  intervals.  Poor  performance  occurred  more 
•frequently  in  the  months  of  December  to  February  and  June  to  August. 

An  additional  approach  of  an  occurrence  frequency  method  is  used  for  determining  the 
performance  of  each  method.  For  each  observation  the  methods  are  ranked  in  descending  order  with 
an  increase  in  prediction  error.  When  the  errors  are  equal  the  respective  methods  are  ranked  equal. 
The  cumulative  count  of  the  respective  methods  is  used  for  determining  the  performance  of  these 
methods.  Also,  the  frequency  is  counted  for  each  method  only  when  the  method  is  the  best  (produces 
the  least  error  among  these  methods)  and  also  for  the  case  when  the  method  is  the  worst  (produces 
the  largest  error).  The  results  for  the  best  and  the  worst  method  approach  (discussed  above)  are 
shown  in  Figure  8.  The  counts  of  the  modified  slope  method  over  the  respective  time  intervals  are 
combined.  The  relative  count  for  performance  for  the  respective  methods  show  that  the  lONCAP, 
the  observed  median  and  the  prior  n  days'  average  methods  are  worse  (short  dotted  lines)  than  the 
averaged  slope  and  the  modified  slope  methods.  The  best  predictions  are  more  frequent  from  the 
modified  slope  and  the  averaged  slope  methods  than  from  the  remaining  methods  (long  dotted  lines). 
The  overall  performance  (continuous  line)  of  the  slope  methods  is  better  than  that  from  the  remaining 
methods.  Considering  all  the  three  groups;  general,  best,  and  worst  performance,  the  modified  slope 
method  with  two  time  intervals  combined,  04-16  for  positive  slope  and  16-04  for  negative  slope,  is 
the  best  method  for  obtaining  f^Fj  predictions  from  the  algorithm  in  Eq.  (1).  The  lONCAP  and  the 
observed  median  (which  is  not  available  in  reality)  methods  are  unable  to  provide  the  f(,F2  predictions 
that  will  satisfy  the  needs  of  the  real  time  operating  systems. 
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Comparison  of  foF2  Prediction  Methods 
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Figure  8.  Performance  of  the  schemes  used  for  the  foF2  predictions. 


4.  CONCLUSION 

The  algorithm  from  Eq.  (1)  proposed  for  improving  the  short  term  real  time  foF2  prediction 
is  better  than  the  other  methods  like  the  lONCAP  and  the  'prior  n  days'  averaged  f^Fj.  Depending 
on  the  circumstances  one  may  truncate  the  algorithm  at  various  terms  on  the  RHS  of  Eq.  (1).  The 
best  results  are  obtained  by  dividing  the  24  hour  period  into  two  intervals  and  using  the  slope 
correction  with  the  appropriate  sign  for  each  interval  to  assure  the  use  of  the  converging  trend  shown 
in  Figure  1.  The  errors  in  the  f^Fj  predictions  are  reduced  from  0.9  MUz  for  the  lONCAP  to  0.4 
MHz  at  high  solar  activity  (1969)  and  from  0.5  MHz  to  0.3  MHz  at  low  solar  activity  (1976).  Thus 
the  proposed  algorithm  improves  the  f^Fj  predictions  by  50  percent  at  high  solar  activity  and  by  40 
percent  at  low  solar  activity.  The  proposed  method,  though  not  perfect,  provides  a  minimum 
improvement  of  25  percent  in  predicting  f^Fj  ^0  percent  of  the  time.  Also,  (Figure  8)  this  method 
has  better  overall  performance,  and  is  the  least  often  the  worst  prediction  method. 

Among  the  various  schemes  considered  here,  the  lONCAP  is  unable  to  provide  a  reasonable 
^2  prediction  for  real  time  operational  systems.  Even  if  the  lONCAP  succeeds  in  reproducing  the 
observed  ^2  median  values,  the  corresponding  improvement  in  the  f^Fj  prediction  is  only  10  percent. 
Prior  'n  days'  averaged  f^Fj  improves  the  prediction  by  14  percent  over  that  from  the  lONCAP 
method,  not  a  very  significant  improvement. 

5.  RELEVANCE  TO  OTH  OPERATION 

Dandekar  and  Buchau^(1986)  have  discussed  the  problem  in  detail  for  the  sunrise  transition 
period.  The  'M'  factor  is  a  ratio  of  the  radar  operation  frequency  to  f^Fj  at  the  midpoint  of  the  radar 
range.  The  radar  range  is  also  a  function  of  the  'M'  factor  and  varies  with  other  F2  layer  parameters; 
Iv  -  peak  height,  ho  -  base  height,  and  y„  -half  width  of  the  layer.  Keeping  all  layer  parameters 
(except  and  radar  frequency  fixed,  the  radar  range  increases  (/decreases)  with  increase 
(/decrease)  in  foF2.  Using  the  equation  relating  range  to  M  developed  by  Appleton  and  Beynon* 
(1940),  Dandekar  and  Buchau*  showed  that  for  the  radar  range  2500-3000  km,  a  2.3  percent  change 
in  M  moves  the  barrier  100  km.  This  corresponds  to  a  6  percent  change  in  M  for  a  sliding  of  the 
barrier  of  500  km,  that  is,  by  half  of  its  width.  In  Figures  5A-E,  changes  in  foF2  larger  than  20 
percent  are  seen  during  the  sunrise  and  sunset  transition  periods.  Dandekar  and  Buchau*  (1986) 
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observed  changes  as  high  as  60  percent  in  fJF^  at  Ottawa  and  Arguello  in  the  winter  season  in  1969. 
If  the  radar  operation  frequency  is  to  be  adjusted  when  the  barrier  slides  by  half  of  its  width,  the 
frequency  has  to  be  changed  every  12  minutes  to  cover  the  f^Fj  change  of  30  percent  and  every  6 
minutes  for  the  change  of  60  percent.  Because  such  rapid  adjustments  of  the  radar  frequency  are 
needed  during  sunrise/sunset  transition  periods,  improved  predictions  of  foFj  with  error  reduced  by 
50  percent  is  very  helpful  to  the  radar  operation. 
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